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ABSTRACT 

We carry out the self-similar solutions of viscous-resistive accretion flows around a 
magnetized compact object. We consider an axisymmetric, rotating, isotheral steady 
I . accretion flow which contains a poloidal magnetic field of the central star. The dom- 

inant mechanism of energy dissipation is assumed to be the turbulence viscosity and 
magnetic diffusivity due to magnetic field of the central star. We explore the effect of 
viscosity on a rotating disk in the presence of constant magnetic diffusivity. We show 
that the dynamical quantities of ADAFs are sensitive to the advection and viscosity 
parameters. Increase of the a coefficient in the a-prescription model decreases the ra- 
dial velocity and increases the density of the flow. It also affects the poloidal magnetic 
k*" | field considerably. 

in 

I ■ 1 INTRODUCTION 

C""^ ' Accretion onto black holes has been intensely studied for the last three decades (see Kato et al. 1998 for a review), and 
several types of models were proposed. Standard accretion disk model of Shakura & Sunyaev (1973) has been very useful 
in interpretation of observations in binary systems and active galactic nuclei, AGN (Pringle 1981), is based on a number of 
^■•^ ' simplifying assumptions. In particular, the flow is assumed to be geometrically thin and with a Keplerian angular velocity 
^ , distribution. This assumption allows gradient terms in the differential equations describing the flow to be neglected, reducing 
K~Ji ' them to a set of algebraic equations and thereby fixes the angular momentum distribution of the flow. For low accretion rates, 
r> '. M, this assumption is generally considered to be reasonable. Since the end of seventies, however, it has been realized that for 
high accretion rates, advection of energy with the flow can crucially modifies the properties of the innermost parts of accretion 
disks around black holes. 

The natural improvement over the Shakura- Sunyeav model was to consider the case when cooling is less efficient than 
viscous heating. This may happen in two cases: either when the disk is extremely optically thick and the radiation is trapped 
for a timescale longer than the accretion timescale( see Abramowicz et al. 1988) or when it is extremely optically thin, a 
regime in which cooling processes are inefficient ( Narayan & Yi 1995b). 

The study of advection accretion flows around low-luminosity black hole candidates and neutron stars are currently a 
very active field of research, both theoretically and observationally (see Narayan et al. 1998 for a review). Observational 
evidences for the existence of low luminosity black holes at the center of galaxies and in the active galactic nuclei AGN 
(Cherepashchuk 1996; Ho 1999) make necessary to revise theoretical models of the accretion disks. Thereby development of 
the subject of advection-dominated accretion flow (ADAF) in recent years lead to global solutions of advection accretion 
disks around accreting black hole systems and neutron stars. In this case, viscously generated internal energy is not radiated 
away efficiently as the gas falls into the potential well of the central mass (as in the standard thin disks model, Shakura & 
Sunyaev 1973) but retained within the accreting gas and advected radially inward (Narayan & Yi 1994, hereafter NY1994) 
and might eventually be lost into the central object or in contrast, a considerable portion of it might give rise to wind onto 
black holes and neutron stars (Blandford & Begelman 1999). By definition, Advection-Dominated Accretion Disks, ADAFs, 
have very low radiative efficiency as a consequence they can be considerably hotter than the gas flow in the standard thin disk 
models (Narayan & Yi 1995 a,b) and therefore they are ultra-dim for their accretion rates (Phinney 1981; Rees et al. 1982). 
On the other hand, since all the internal energy is stored as thermal energy, the gas becomes extremely hot and the kinetic 
temperature of ions approaches the virial limit what means that even at high initial angular momentum, the disk becomes 
very thick, forming practically a quasi-spherical accretion flow (Narayan & Yi 1995a). A general description of an advective 
accretion flow around a compact star was put forward by Narayan & Yi (1994, 1995a): they parameterized the degree of 
advection with one parameter /, defined as the ratio between the thermal energy stored in the disk and advected the central 
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object ( not radiated), and the total thermal energy generated by viscosity at each radius. The general result obtained both 
from self-similar solutions and numerical calculations is that high advection (/ ~ 1) produces a hotter, thicker disk with a 
larger infall radial velocity and conversely a sub-Keplerian circular velocity. 

A remarkable problem arises when the accretion disk was threaded by magnetic field. There are good reasons for believing 
that the magnetic fields are important in accretion processes in ADAFs. Some authors tried to study magnetized accretion 
flows analytically. For example ,Kaburaki (2000) presented a set of analytical solutions for a fully advective accretion flow in 
a global magnetic field. Shadmehri (2004), hereafter SH2004, extended this analysis by considering a non-constant resistivity. 
He obtained a set of self-similar solutions in spherical coordinates that described quasi-spherical magnetized accretion flow. 
The full account of the processes, connected with a presence of magnetic field in the flow, is changing considerably the picture 
of the accretion flow. 

In ADAF models, energy dissipation in the accretion flow can be assumed to be due to turbulent viscosity and electrical 
resistivity. Under some conditions, it is important that we consider the effect of resistivity on accretion flows. Kuwabara et al. 
(2000) showed the results of global MHD simulations of an accretion flow initially threaded by large-scale poloidal magnetic 
fields including the effects of magnetic turbulent diffusivity. They found the importance of strength of magnetic diffusivity 
when they studied it in magnetically driven mass accretion. They pointed out that the mass outflow depends on the strength 
of magnetic diffusivity, so that for a highly diffusive disk, no outflow takes place. Thereby in this paper, we want to explore 
how the structure of a steady state thick disk depends on its resistivity and viscosity so we pursue SH2004's work, resistive 
disks, when accretion flow experiences the rotation as well as viscosity dissipation. We consider ADAFs with the pure inflow 
and investigate the effect of viscosity on some physical quantities of the flows such as the radial and angular velocities, the 
density and the magnetic field flux. In order to study the dynamics of these flows, several simplified assumptions must be 
made in the analysis. The fluid is treated at least approximately as non-relativistic and also a poloidal model is adopted for 
electromagnetic field in which, it has a poloidal component in the disk. However, we will present self-similar solutions for 
viscous-resistive ADAFs. 

This paper is organized as follows. Section 2, we present the equations of magnetohydrodynamics as the basic equations. 
General principles are presented in section 3. We show that the equations can be solved using the self-similar method and the 
numerical solutions are discussed in section 4 followed by results in section 5, summery and conclusion in section 6. 



2 THE BASIC EQUATIONS 

As we stated in introduction, we are interested in constructing a model for describing accretion disks in global magnetic fields. 
The macroscopic behavior of such flows can be studied by MHD equations. For simplicity the self-gravity of the disks and the 
effect of general relativity have been neglected. The flow is described in terms of the flow-frame time derivative or co-moving 
derivative, i.e. that defined as: = + (u • V). So, we can describe the accretion flows by the fundamental governing 
equations which are written by the equation of continuity: 

^+pV-u = 0, (I) 
the equation of motion: 

P^- = -Vp - pV$ + fiV 2 u + {fit + |) V(V • u) + i- J x B, (2) 
the equation of energy: 

Q+-Q- = Q adv , (3) 
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and the Maxwell's equations: 

V ■ B = 0, (4) 
— = V x (u x B) + n V 2 B, (5) 

where p is the density of the gas, p the gas pressure, e the internal energy, u the flow velocity, B the magnetic field, J = VxB 
the current density, r\ the magnetic diffusivity in which for simplicity it is assumed to be a constant parameter (see, e.g., 
Kaburaki 2000), /i and /lb are the shear and bulk viscosities, Q adv represents the advective transport of energy and is defined 
as the difference between the viscous heating rate Q + and radiative cooling rate Q~ . We neglect self-gravity so that $ is 
assumed to be due to a central object. Also we neglect radiation pressure in the equations because in the optically thin 
ADAFs, P 9a3 > P rad . 

We employ the parameter / = (e.g., NYf994) to measure the degree to which the accretion flow is advection- 
dominated. In general, it will vary with r and depends on the details of the heating and cooling processes. For simplicity, it 
is assumed a constant. So, for advection-dominated flows we have Q adv ~ Q + >• Q~ . This corresponds to an optically thin 
ADAF where the viscous energy is stored in the gas as internal energy and the amount of cooling is negligible compared 
to the heating. In this case, the accreting gas has a very low density and accretion rates are low, sub-Eddington (Ichimaru 
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1977; Rees et al. 1982; NY1994; Abramowics et al. 1995). Correspondingly, for radiative-cooling-dominated flows we have 
Q + ~ Q~ 3> Q adv . This corresponds to a cooling-dominated flow where the viscous energy is released in the gas as radiative 
energy and the amount of energy advected is negligible. The optically thick Shakura-Sunyaev disks correspond to this case. 
Now, we formulate the basic equations (HJ-© in spherical polar coordinates as follows: 

9 P , 1 d I 2 1 , 1 d I \ n !a\ 

di + ^Tr^ PUr) + - W (pue)=Q, (6) 
The three components of the momentum equations give (e.g., Mihalas & Mihalas 1984): 
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also, the r-component of the induction equation 
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3 GENERAL PRINCIPLES 

As stated, we used a total of eight partial differential equations governing the non-self gravitating flow. These equations relate 
15 dependent variables: p, p, e, (i, fib,ri and the components of u, J and B. The equations must be closed by specifying suitable 
prescriptions for the viscosity and for the thermodynamics. Thus for the set of equations adopted, we make the following 
standard assumptions: we consider a steady (J^ = 0), rotating, axisymmetric (^ = 0), quasi spherical accreting flow (so it 

is convenient that we use spherical coordinates (r, Q,<f>) in our discussion) with Keplerian angular velocity fifc(r) = (-^i^) 5 
around a central object and with a purely poloidal magnetic field threading the disk. We assume that the fluid can be treated 
at least approximately as non-relativistic. 

The kinematic viscosity coefficient, v = — , is generally parameterized using the a-prescription (Shakura-Sunyaev 1973), 

v = ac s H, (14) 

where H = -Sf- is known as the vertical scale height , c s = .f^ is the isothermal sound speed and the dimensionless coefficient 

* L k VP 

a is assumed to be independent of r. Also it is important that we consider the effect of resistivity on accretion flows. So, 
we introduce the parameter r\ as the magnetic diffusivity and insert it as a constant parameter in our equations. Both the 
kinematic viscosity coefficient v and the magnetic diffusivity r\ have the same units and are assumed to be due to turbulence 
in the accretion flow. Thus it is physically reasonable to express r\ such as v via the a-prescription of Shakura-Sunyaev (1973) 
as follows (Bisnovatyi-Kogan & Ruzmaikin 1976), 

r) = r] c s H. (15) 

where the dimensionless coefficient rjo is assumed independent of r. Bisnovatyi-Kogan and Ruzmaikin (1976) introduced rj 
and proposed that rj ~ a. Substituting H and c s into the relations (|14[) and (|15|) . we find them locally proportional to the 
pressure: 

" = (16) 

V = Vo-^-. (17) 
pih 

Their ratio is one definition of the magnetic Prandtl number as follows, 

P m = V - = ^A. (18) 
«7 V 

where ut is the random velocity of diffusing particle and I is its mean free path. For a fully ionized hydrogen plasma Prandtl 
number is very large compared with unity due to the large length scale I of the disk. In the following we consider conditions 
where Pm ^ 1. That means, viscous and resistive forces can be contributed in the energy dissipation similarly when it is 
equal to unity and resistive forces are dominate when P m is smaller than unity. 

To determine thermodynamical properties of the flow in the energy equation (|10[) . we require a constitutive relation as a 
function of two state variables. Therefore we choose an equation for the internal energy as e — p ^_^ where V is the ratio of 
specific heats of the gas. 

To satisfy V ■ B = 0, we may introduce a convenient functional form for the magnetic field. Owning to the axisymmetry, 
the magnetic field can be written as 
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B = B p (r, 9) + B<p(r, #)e^ (19) 

The effect of magnetic diffusivity on magnetically driven mass accretion was studied by Kuwabara et al. (2000). They 
showed that the effects of resistivity are that magnetic field lines do not rotate with the same angular speed as the disk matter 
and thus it suppresses the injection of magnetic helicity and magneto-centrifugal acceleration. So, neglecting the toroidal 
component of field, B^, we can express the poloidal component, B p , in terms of a magnetic flux function ^(r, 9): 

B = B p (r, 9)4vx (-^ e «) . (2°) 
2tv \r sm 9 J 

in which satisfies V • B = V ■ B p = 0. Magnetic flux function ^(r, 9) is related to the magnetic vector potential by 9 — rA^ 
with A$ the toroidal component of the vector potential. The magnetic flux contained inside the circle r=constant, 9 =constant 
is, 

r 

B g (r',6)2irr'dr' = 2-K^{r,9) (+constant), (21) 
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Since B p ■ = 0, $ labels magnetic lines or their surfaces of revolution, magnetic surfaces. 
Similarity, one can write the flow velocity in the form 

u = u p (r, 9) + (r, 0)e Kj> , (22) 

subsequently, we show that the poloidal component of velocity has only a radial component u p = u r . We take u r to be a 
negative, since we want to consider infall of material. 

It is clear that the basic equations are nonlinear and we cannot solve them analytically. Therefore, it is useful to have a 
simple means to investigate the properties of solutions. This is most easily done in terms of a set of dimensionless parameters 
which can be expected to be similar at all times. Here, one can employ the method of self-similar to fluid equations. 



4 SELF-SIMILAR SOLUTIONS 

To better understanding the the physical processes of our viscous-resistive ADAF accretion disks, we seek self-similar solutions 
of the above equations. The self-similar method is familiar from its wide applications to the full set of MHD equations. As 
long as we are not interested in boundaries of the problem, such solutions can accurately describe the behavior of the solutions 
in an intermediate region far from the radial boundaries. 

Writing the equations in non-dimensional forms, that is, scaling all the physical variables by their typical values, brings 
out the non-dimensional variables. We can simply show that a solution of the following forms, satisfy the equations of our 
model: 

p(r,9) = Po p(9)(r/r o y 3/2 , (23) 
p(r,9)=p P(9)(r/r )- 5/2 , (24) 
u T {r,9)=rQ, K (r)U{9), (25) 
u 9 {r,9) = rsmen K (r)n(9), (26) 

*<'••>- -S'^""- < 28 » 

where p , Po, B a and r a provide a non-dimensional form for our equations. Also ug = is considered, since the most 
important assumption we made is that the flow is steady (Jj = 0) and also with respect to above solutions we find r 2 pu r is 
independent of r but it is a function of 9). It represents mass accretion rate in a given 9. If we integrate this over the angle 9, 
we obtain the net mass accretion rate 



M = -2hL J pU sin 9d9 (29) 
where M a = 2-np rl s [^-. 

We adopt Q-prescription (|16p so that p — vp — a jf-. The bulk viscosity is not usually discussed in the context of accretion 
flows. Thus we assume pb = 0. Substituting the above solutions in the equations (f7ll- (|13[) . we obtain a set of coupled ordinary 
differential equations in terms of 9 for the three components of the equation of motion: 
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the energy equation becomes: 
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The above equations constitute a set of ordinary differential equations for functions U(9), 9(9), P(0) and 0,(6) as follows: 
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constant. Also we can obtain p as follows, 
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b tj^- Integrating the last equation, and doing some simplifications we have f2<E' = K, where K is an arbitrary 
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From equation (|40p we can find: 
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Comparing equations (|37|l and (|42|) we can eliminate $ and then the result can be compared with equation (|38[) , finally we 
have: 
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Equation (|43|l is an ordinary first order differential equation for = 48 , and it has two roots: 
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substituting the above solution in the main equations of the system, (|36|I - (|37| I- (I39I) . and eliminating $ we obtain: 
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Equations (J44J)- (J4TJ) constitute a system of ordinary non-linear differential equations for the four self-similar variables, 
Q, P, U, p. Indeed, the behavior of the solutions depends on boundary conditions which are supposed based on some physical 
assumptions such as symmetry with respect to the equatorial plane, etc. 

There are many techniques for solving these nonlinear equations. Analytical methods can yield solutions for some simplified 
problems. But, in general this approach is too restrictive and we have to use the numerical methods. Here, one can employ 
the method of relaxation to the fluid equations (Press et al. 1992). In this method we replace ordinary differential equations 
by approximate finite-difference equations on a grid of points that spans the domain of interest. The relaxation method 
determines the solution by starting with a guess and improving it, iteratively. Based on it, this system of equations can be 
solved for all unknowns as a function of 9, once if we have a complete set of boundary conditions which put some physical 
constraints on the flow. The boundary conditions are distributed between the equatorial plane, 9 = ? and the rotation axis, 
9 = 0. The boundary conditions at 9 — are: 

du dn dp TT 

■d9 = d9 = M = M=° ' U = °> (48) 
and in this method the boundary conditions at 9 — ^ are: 

dU _ dQ _ dP _ 
lW~liB~~dB~lW~ •' 



(49) 



We have found condition %z 



by the relation QA>^ — K for both two boundaries which demand the magnetic flux 
enclose by the polar axis goes to zero. So field line thread the equator vertically. We have a limitation to reach 9 = 
numerically. So we try to use boundary condition very close to polar axis with very small value of 9. We can obtain physical 
solutions if we consider the big value of Q and so VP ~ at 8 — e. Also, at 9 = ? the equations (|44[1 and (|47[) gives: 
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The boundary conditions on the above equations require that variables are assumed to be regular at the endpoints. Also the 
net mass accretion rate (|29l) provides one boundary condition for p: 



p{9) U{9) sin ( 
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For solving the MHD equations numerically, we need these boundary conditions. Using these boundary conditions require to 
choose K, e ,5,7 properly. So we are going to have a astrophysical approximations for this variables. 

For our illustrative parameters, we consider approximate equipartition of magnetic and kinetic energies for matter flowing 
into a neutron star. In quasi-spherical accretion flows, approximate equipartition proposed by Shwartsman (1971), 

B 2 1 

8^ ~ 2 pU " (50) 

where u r corresponds to an Alfven speed ua = yf= — Uff (free-fall speed). The density of the infalling materials can be 
reasonably approximated by free-fall behavior. In the accretion disks around the magnetized compact object, due to the 
complicated nonlinear effects from the interaction between accretion materials and the magnetic fields of the compact object, 
the mechanisms that led to the magnetic field chandeled accretion flow still remain unsolved. However, it is widely accepted 
that due to strong gravity of the compact object, a very steep and supersonic flow hits the magnetosphere boundary and 
deforms its structure drastically. The infalling flow is halted at distance r ~ ta where the magnetic energy density balances 
the flow kinetic energy, or ua — us- Here ua is calculated from free fall density and velocity: 



M 2GM 

PS = ~. 5 , US = \ (51) 

where M is the mass accretion rate and ta Alfven radius. Interestingly, the numerical simulations (Paatz & Camenzind 1996, 
Koldoba, Lovlace & Ustyugova 2002) show that at this distance r ~ ta the material funnels through the magnetosphere and 
fall onto the magnetic poles. The funnel flow is found to reach the star's magnetic poles with the velocity close to that of 
free-fall velocity. 

The magnetic configuration outside of a compact star, such a neutron star, is a dipole field then in spherical polar 
coordinates we have 
t, ,-, r>3 ( 2 cos sm9 \ 

B = B v R + — 66 j (52) 

where R is the radius of the neutron star and B p is the magnetic field at the magnetic poles. Thus in the equatorial plane 
(0 — j), the equation (|52[) becomes B ~ B " s . So we obtain, 

4nr 2 A u ff 8nGMr%' y ' 

The equation (153[) implies that the infalling mass close to the equatorial plane of a neutron star is stopped at a radial distance, 
r A = B^R^M-^^GM)- 1 ' 7 (54) 

for M ~ M Q , M ~ 1.85 x 10 16 ffs _1 (M < Mem = 1.39 x lO 18 ^^ -1 ), B p ~ 10 12 G and R ~ 10 6 cm, the equation (54j) gives 

r ~ ta ~ 5.2 x 10 cm. Thus we can obtain B = B p ~ 10 6 G and p a ~ 3.6 x 10 _10 5 cm . On the other hand, the gas 

pressure, density and temperature are related by 

Po = (55) 

p, a m H 

where p, a is the mean molecular weight, k Boltzman constant and mn ~ m p is the mass of the hydrogen atom. We assume 
the accreting materials to be mainly hydrogen and to be nearly ionized. For accreted materials, the potential energy released 
is GM ; M °" and the thermal energy is 2 x |fcT (T vir = 2k.); therefore, T th = GA 3*" acc ~ 5.5 x W 9 K we obtain p a ~ 



3kr 

3.2 x 10 8 dyne cm~ 2 . Consequently, in turn it can be reduced that j3 = 0.01, C\ = 0.8 and C2 = 2 x 10 3 . 

Imposing two boundary conditions (|48[) and (J49J to the equations (|44I) - (I47I ). we obtain numerical solutions for the flows 
with a fixed value of r) = 0.1, T = 4/3 and a = 0.01, 0.05 and 0.1 with a sequence of the increasing values of parameter f or 
decreasing cooling. Figures 1-5 show our results. 



5 RESULTS: TYPICAL SOLUTIONS 

We have obtained numerical solutions of equations (|44p ~(|47 p for variety of values of the viscosity parameter a and advection 
parameter /. Figures 1-4 show a typical sequence of solutions correspond to a — 0.01,0.05,0.1 and / = 0.1,0.3,0.5,0.7. 
The solutions may be considered either as flows with a fixed value of viscosity parameter and with a sequence of increasing 
advection parameter or decreasing cooling (Figures 1 and 3) or a sequence of different values of viscosity parameter a in a 
fixed advection regime (Figures 2 and 4). 

The six panels in Figure 1 show the variations with respect to the polar angle 9 of various dynamical quantities in the 
solutions. The top left panel displays the dimensionless radial velocity U(8) as a function of 8. The velocity is zero at 9 — 
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Figure 1. The self-similar solutions of radial velocity U (left panels) and density p (right panels) as a function of polar angel 9 
corresponding to for Top: a=0.01, Middle: «=0.05 and Bottom: a=0.1. 

(this is a boundary condition) and maximum at = ^. So the maximum accretion velocity is at equatorial region and on 
the polar axis there is no mass inflow. As we expected the velocity is sub-Keplerian. We find that in the boundary, U(8) is 
essentially independent of advection parameter /. But in the intermediate, the radial velocity is modified by /; in the SH2004 
solutions, two distinct regions in the U(9) profile could be recognized. The bulk of accretion occurs from equatorial plane at 
9 = | to a surface at 6 = 9 a , inside of which the radial velocity is zero. While NY1994 solutions there is no zero inflow in 
< 6 < . Our solutions show that in a given 6 the radial velocity is increased when we increase the advection parameter. The 
middle and bottom left panels display the radial velocity in a = 0.05, 0.1 for a sequence of advection parameters respectively. 

The top right panel shows profile of the density p{9). The density contrast in the equatorial and polar regions increase 

with decreasing advection parameter /. The density grows and becomes concentrated toward the equatorial plane. For a given 

a, solutions with small values of / behave like standard thin disks, as might be expected since these solutions correspond 

to / — ■» and so advect very little energy. In the opposite advection-dominated limit, which corresponds to / — » 1, our 

solutions describe nearly spherical flows which rotate at much below the Keplerian velocity. This is demonstrated in Figure 

5 where we display iso-density contours in the meridional plane. The middle and bottom right panels display the density 

profiles in q = 0.05,0.1 for a sequence of advection parameters respectively. This advection-dominated solutions have very 

similar properties to the approximated solutions derived by NY1994 and SH2004. The results show that these quantities 

are very sensitive to the advection parameter. For low advection regimes, / = 0.1,0.3, we have to separate the regime for 

radial inflow. But for high advection regimes / = 0.5, 0.7, the radial velocity is non zero around the pole. The behavior of 

the solutions, self-similar radial velocity and density profile, are demonstrated in Figure 2 for different values of advection 

parameters with variation of viscosity parameters. The panels in Figure 2 show that For a fixed /, the density maximum 

© 2007 RAS, MNRAS 000, Q3-?? 
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Figure 2. The self-similar solutions of radial velocity U (left panels) and density p (right panels) as a function of polar angel 9 
corresponding to T = 4/3, »?o = 0.1, /3 = 0.01 and a = 0.01, 0.05, 0.1 for / = 0.1, 0.3, 0.5, 0.7 from Top to Bottom, respectively. 



increases with increasing viscosity parameter while the radial velocity decreases with increasing it. In this case, increasing the 
viscosity parameter corresponds to increasing heating mechanisms, so in a fixed advection regime, there are more energy to 
advect into central star. The top Panels in Figure 5 Show that the disk to be thick. The solution with the same / but with 
different values of a are virtually indistinguishable from one another. There are probably, more significantly variations when 
a exceeds above 0.1. Recently King et al. 2007 assart that in a thin and fully ionized disk the best observational evidence 
suggest a typical range a ~ 0.1 — 0.4 where relevant numerical simulations tend to drive estimates for a which are an order 
of magnitude smaller. However, such large values of a are probably unlikely (eg. Narayan, Loab, & Kumar 1994; Hawley, 
Gammie, & Balbus 1994), and so we have not explored this region of the a parameter. The U{6) and p(9) profiles both peak 
at 6 = Therefore, in all our solutions the bulk of accretion occurs along the equatorial plane and accretion rate goes to 
zero along the rotational pole. An interesting feature of these solutions, as already mentioned, for low advection ,/, we have 
a thin disk; in all cases there is a low density with a higher temperature corona above the disk. 

Figure 3 displays the magnetic flux function and the angular velocity for different values of advection parameter in a fixed 
viscosity parameter. Integration of equation (|40|l yields f^ 2 = K, where K is an arbitrary constant. So angular velocity and 
magnetic flux function should have opposite behavior along 9 direction. We plot them with choosing proper input parameters 
which introduced in the end of the last section. Figure 3 shows that magnetic flux function, ^{0), varies by only ~ 50 
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Figure 3. The self-similar solutions of magnetic flux (left panels) and angular momentum Q (right panels) as a function of polar angel 
9 corresponding to T = 4/3, r) = 0.1, /3 = 0.01, / = 0.1,0.3,0.5,0.7 for Top: o=0.01, Middle: a=0.05 and Bottom: a=0.1 in rotating 
accretion flows. 



percent. The magnetic flux function increases by increasing advection in accretion disk in a fixed viscosity. The behavior of 
angular velocity is exactly opposite. Figure 4 displays the behavior of the magnetic flux and angular velocity for different 
values of viscosity parameter in a fixed advection. The solutions implies that in a fixed, low /, the effect of different a is 
indistinguishable. 



6 SUMMERY AND CONCLUSION 

The main aim of this investigation was to obtain axisymmetric self-similar advection-dominated solution for viscose-resistive 
accretion flow. We have presented the results of self-similar solutions of the effect of the viscosity and rotation on magnetically 
driven accretion flows from a flow threaded by poloidal magnetic fields where only serious approximation we have made is 
the use of an isotropic a viscosity and constant magnetic diffusivity. Attention has restricted to flow accretions in which self- 
gravitation is negligible. We included the magnetic diffusivity so that its value was constant throughout our analysis. Using the 
basic equations of fluid dynamics in spherical polar coordinates (r, 6, </>), we have employed the method of self-similar for thick 
discs to derive a set of coupled differential equations which govern the dynamics of the system. We then solved the equations 
by the method of relaxation by considering boundary conditions and using a-prescription (Shakura & Sunyaev, 1973) in order 
© 2007 RAS, MNRAS 000, []>?? 
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Figure 4. The self-similar solutions of magnetic flux ^ and angular momentum Q as a function of polar angel 8 corresponding to 
T = 4/3, rjo = 0.1, A) = 0.01 and a = 0.01, 0.05, 0.1 for / = 0.1, 0.3, 0.5, 0.7 from Top to Bottom, respectively. 



to extract some of the similarity functions in terms of the polar angle 8. Figures are considered for a = 0.01, 0.05, and 0.1 so 
that for any a we considered / = 0.1, 0.3, 0.5 and 0.7. 

We showed that the radial and rotational velocities are well below the Keplerian velocity. The Bulk of accretion with 
nearly constant velocity occur in the region which extend from equatorial plane to a given 9 which highly depends on advection 
parameter /. In a non-advective regime, low /, we have a standard thin accretion disk but for high / the accretion is nearly 
spherical. The geometrical shape of the flow is determined by the amount of viscosity and advection, in a fixed magnetic 
diffusivity. The accretion disk with efficient cooling (/ — > 0) has low-density, high temperature corona which implies that the 
regular thin disks may be accompanied by advection-dominated corona which can drive low-density wind. 

Our results show the flow of ionized accretion materials is not disklike in morphology. The closest along our solution in 
the accretion literature is Bondi (1952) spherical accretion. Our flows differ in important ways from Bondi problem. The gas 
in our model rotates and has viscose interaction through which angular momentum is transported outward. Also, our flow has 
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Figure 5. Isodensity contours in the meridional plane for T = 4/3, rjo = 0.1, /3o = 0.01 for Top: /=0.5 and a=0.01, 0.05, 0.1, Bottom: 
a=0.01 and /=0.1, 0.3, 0.5, 0.7. 

magnetic interaction with dipole magnetic field of the central star which it can redistribute the angular momentum within 
the accretion disk. The angular velocity is significantly sub-Keplerian, and this may have an important role in spin-up of 
accreting stars. Stars which their spins grow by interacting with accretion materials are likely to reach a steady state with 
a rotation rate below the break-up limit. This advection-dominated solution may be a good solution for angular momentum 
problem in the star formation. 

However, our results improve the physics of advection-dominated accretion disks. It is important that the magnetic 
diffusivity can modify the dynamical quantities of the disks. We developed NY1994 solutions to a realistic model for ADAFs 
by adding magnetic diffusivity. But in future we are going to investigate the effect of non-constant magnetic diffusivity, rio- 
Several developments can be investigated to reach a much realistic description of the physics of accretion disks around the 
magnetized compact objects. 
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This document presents a number of hints about how to set up your Science 
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paper in Mj<X . We provide a template file, scif ile . tex, that you can use 
to set up the ETgX source for your article. An example of the style is the special 
^ ■ {sciabstract} environment used to set up the abstract you see here. 

o Introduction 

o 



In this file, we present some tips and sample mark-up to assure your ETgX file of the smoothest 
possible journey from review manuscript to published Science paper. We focus here particularly 
on issues related to style files, citation, and math, tables, and figures, as those tend to be the 
biggest sticking points. Please use the source file for this document, scifile.tex, as a 
template for your manuscript, cutting and pasting your content into the file at the appropriate 
places. 

Science's publication workflow relies on Microsoft Word97. To translate ETgX files into 
Word97, we use an intermediate MS-DOS routine (?) that converts the T^X source into HTML. 
The routine is generally robust, but it works best if the source document is clean ETgX without 

1 



a significant freight of local macros or . sty files. Use of the source file scif ile . tex as a 
template, and calling only the .sty and . bst files specifically mentioned here, will generate 
a manuscript that should be eminently reviewable, and yet will allow your paper to proceed 
quickly into our production flow upon acceptance (?). 

Formatting Citations 

Citations can be handled in one of three ways. The most straightforward (albeit labor-intensive) 
would be to hardwire your citations into your ETgX source, as you would if you were using an 
ordinary word processor. Thus, your code might look something like this: 

However, this record of the solar nebula may have been 
partly erased by the complex history of the meteorite 
parent bodies, which includes collision-induced shock, 
thermal metamorphism, and aqueous alteration 
({\it 1, 2, 5 — 7\/}) . 

Compiled, the last two lines of the code above, of course, would give notecalls in Science style: 

. . . thermal metamorphism, and aqueous alteration (1, 2, 5-7). 

Under the same logic, the author could set up his or her reference list as a simple enumera- 
tion, 

{ \bf References and Notes} 
\begin { enumerate } 

\item G. Gamow, {\it The Constitution of Atomic Nuclei 
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